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Abstract: Investigating the relationships between vegetation dynamic and edaphic factors provide 
management insights into factors affecting the growth and establishment of plant species and vegetation 
communities in saline areas. The aim of this study was to assess the spatial variability of various 
vegetation communities in relation to edaphic factors in the Great Salt Desert, central Iran. Fifteen 
vegetation communities were identified using the physiognomy-floristic method. Coverage and density of 
vegetation communities were determined using the transect plot method. Forty soil samples were collected 
from major horizons of fifteen profiles in vegetation communities, and analyzed in terms of following soil 
physical and chemical characteristics: soil texture, soluble Na* concentration, sodium adsorption ratio 
(SAR), electrical conductivity (EC), pH, organic matter content, soluble Mg?* and Ca?* concentrations, 
carbonate and gypsum contents, and spontaneously- and mechanically-dispersible clay contents. 
Redundancy analysis was used to investigate the relationships between vegetation dynamic and edaphic 
factors. The generalized linear method (GLM) was used to find the plant species response curves against 
edaphic factors. Results showed that plant species responded differently to edaphic factors, in which 
soluble sodium concentration, EC, SAR, gypsum content and soil texture were identified as the most 
discriminative edaphic factors. The studied plant species were also found to have different ecological 
requirements and tolerance to edaphic factors, in which Tamarix apbylla and Halocnemum strobilaceum were 
identified as the most salt-resistant species in the region. Furthermore, the presence of Artemisia sieberi was 
highly related to soil sand and gypsum contents. The results implied that exploring the plant species 
response curves against edaphic factors can assist managers to lay out more appropriate restoration plans 
in similar arid areas. 


Keywords: desert; generalized linear method; mechanically-dispersible clay; ordination; plant species response cutves; 
spontaneously-dispersible clay 


1 Introduction 


Understanding how vegetation coverage varies in response to edaphic factors, especially in arid 
environments, can serve as a useful guideline for effective land use and identifying appropriate 
species for desert restoration programs (Mokhtari-Asl et al., 2008; Sheikhzadeh et al., 2019). At 


"Corresponding author: Hossein BASHARI (E-mail: hbashari @cc.iut.ac.ir) 
Received 2019-04-15; revised 2020-10-15; accepted 2020-12-23 


© Xinjiang Institute of Ecology and Geography, Chinese Academy of Sciences, Science Press and Springer-Verlag GmbH Germany, 
part of Springer Nature 2021 


http://jal.xjegi.com; www.springer.com/40333 


JOURNAL OF ARID LAND 


large scales, climatic factors exert the greatest influence on the growth and distribution of plant 
species, while soil physical and chemical properties appear to play a more important role at local 
scales (Kent, 2011; Khatibi et al., 2017). Halophytes are salt-tolerant plants that complete their 
growth cycles thoroughly in saline soils (Flowers et al., 1977; Agarwal et al., 2020). The growth 
and reproduction performance of these plant species considerably differ across the salinity 
gradient (Ungar, 1967). Saline soils are characterized by a low plant diversity with a small 
number of (or even a single) plant species, whose distributions are associated with their own 
biological characteristics and the region's edaphic properties (Abbas et al., 2016, 2017). Despite it 
is believed that edaphic properties such as electrical conductivity (EC), acidity and water table 
level substantially influence the vegetation communities of halophytes, more assessments need to 
be conducted to gain a clear understanding of plant distribution in saline soils (Salama et al., 2013; 
Gonzalez-Alcaraz et al., 2014; Abbas et al., 2016). 

Identifying highly resistant species to harsh environments and soil/environmental stresses is 
important for succeeding land management planning and combating desertification practices 
(Gonzalez-Alcaraz et al., 2014). Plant species response curves have substantial implications for 
the interpretation of the relationships between plant species and edaphic conditions. These curves 
can be derived by statistical techniques such as the generalized linear model (GLM) and 
generalized additive model (GAM) as well as data from multivariable analysis. During the past 
decades, many attempts have been made to assess the growth performance of plant species against 
edaphic factors (Magurran, 2004; Park et al., 2020). Most studies in this field have selected small 
geographical extents with a relatively narrow range of habitat parameters that may potentially 
lead to biased interpretation of the species ecological niche (Piernik, 2012), while neglecting to 
determine how significant these findings may be over vast environmental gradients in desert and 
saline regions. 

The flora of arid regions in central Iran belongs to the Irano-Turanian ecoregion and contains a 
large number of endemic plant species threatened by a heightened risk of extinction, mainly due 
to adopting improper policies such as dam construction that diverts water flow and cuts water 
supply to central Iran. Such improper practices contribute to increasing the area of saline soils and 
accelerate desertification in these fragile areas. It is estimated that around 24x10° hm? (15%) of 
land in Iran is covered with salt-affected soils (Piernik, 2012; Allahgholi and Asri, 2014). 
Although the association between soil properties and vegetation coverage was addressed well in 
some arid and saline areas (Abbas et al., 2016; Koull and Chehma, 2016; Salama et al., 2016; Cao 
et al., 2020), it has not yet received the desired attention in some other parts of the world such as 
central Iran, which suffer from limited accessibility and extreme environmental conditions. 

This study aimed to investigate the relationships between vegetation coverage and edaphic 
factors using various multivariate analyses. The knowledge on these relationships and 
extrapolating the results to other similar regions can provide practical solutions to rehabilitating 
and restoring degraded rangelands in arid and desert areas. We also investigated the contribution 
of soil dispersible clay (i.e., an index of aggregate instability), which plays a pivotal role in soil 
erodibility and plant establishment. Another attempt was made to study and classify the plant 
species response curves in relation to edaphic factors in the desert area of central Iran. The results 
of this research will contribute for protecting the endemic species from the risk of extinction and 
preventing of desert encroachment on agricultural lands. 


2 Materials and methods 


2. Study area 


The study area of this research comprises the Hoz-e Soltan Playa, the Moreh and Masileh plains 
in Qom Province, which is located in the northwestern edge of the Great Salt Desert, central Iran 
(34?45'-34555'N, 51?15'-51?30'E; Fig. 1). The study area covers an area of 16x10* hm? and has 
an arid and semi-desert climate based on the De Martonne climate classification method (de 
Martonne, 1926). The coldest and hottest months are January and June in this area with mean 
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temperatures of —1°C and 40°C, respectively. The mean minimum and maximum relative 
humidities are 5% and 98%, respectively. The monthly evaporation is the maximum in June (i.e., 
493 mm). The annual potential evapotranspiration is 1958 mm (IMO, 2015). Precipitation has a 
large spatiotemporal variability in the study area. From south and west to north and east, 
temperature increases and precipitation and the number of freezing and cloudy days decrease. 
With a gentle slope, the plains of the region are extended inward to the salt lakes causing the land 
to be formed mostly by alluvial deposits with varying soil texture. Some parts of the study area 
have high soil infiltrability whose subsurface water resources are near to the soil surface (Asrari 
et al., 2012). 
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Fig. 1 Location of the study area (a) and sampling sites (b) in the Great Salt Desert, central Iran 


2.2 Data collection and analysis 


Vegetation communities were firstly distinguished by overlaying aerial photographs with 
1:50,000 topographical maps. Afterwards, the field-based physiognomy-floristic method was 
employed to better identify vegetation communities. In total, fifteen vegetation communities 
(sampling sites) were identified across the study region and termed as S,—Sjs in this study (Fig. 
1b). A random systematic sampling design was used to collect soil and plant samples from each 
vegetation community. We used the minimal area method to determine the optimum plot size for 
studying vegetation characteristics (Mueller-Dombois and Ellenberg, 1974). Plot sizes varied 
from 4 m? (2 mx2 m) in the Artemisia community to 25 m (5 mx5 m) in the Tamarix 
community. Moreover, we determined the number of plots (N) following the initial sampling 
based on the changes in vegetation coverage using the Cochran's method (Eq. 1) (Cochran, 1977): 


N = (P x s) / (k x x), (1) 
where t is the value of f-student with n—1 degrees of freedom at the probability level of 5%; s is 


the sample variance; xis the mean; and k is the acceptable error value. 

Three 50-m transects, each comprising ten appropriately sized plots were established in each 
vegetation type. Number of plant species and their coverage were quantified within each plot. In 
order to determine the physical and chemical properties of soils, we excavated 1-m vertical soil 
profile in each vegetation community. We collected soil samples from each soil horizon based on 
visual inspection (i.e., 40 soil samples from 15 soil profiles). Samples were air-dried and passed 
through a 2-mm sieve mesh. In the laboratory, soil samples were analyzed to determine soil 
texture by hydrometer method, and pH and electrical conductivity (EC) were determined by pH 
and EC meters, respectively. Soil calcium carbonate equivalent (CCE) was measured by NaOH 
back-titration method, gypsum by acetone method, soluble Ca^ and Mg? concentrations by 
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titration method, soluble Na* concentration by a flame photometer, and organic matter (OM) 
content by Walkley-Black method (Carter and Gregorich, 2008). Soil aggregates instability was 
characterized using spontaneously-dispersible clay (SDC) and mechanically-dispersible clay 
(MDC) contents. SDC characterizes the soil dispersible clay without any external force and is 
important in sodic soils, but MDC accounts for the soil dispersible clay after application of an 
external force, e.g., shaking (Rengasamy et al., 1984; Czyz and Dexter, 2015). SDC and MDC 
were determined by sedimentation (pipette) method. We calculated SDC and MDC as 
percentages of total clay in order to characterize the dispersible contents independent of soil 
texture (total clay). 

Weighted mean properties of soil horizons (Weaver et al., 1994) were used to investigate the 
effect of soil on vegetation coverage. We analyzed the relationships between vegetation coverage 
and edaphic factors using various multivariate analyses based on the canopy cover percentage and 
soil „properties (soil texture, soluble Na*, sodium adsorption ratio (SAR), EC, pH, OM, soluble 

Mg * and Ca^* concentrations, carbonate and gypsum contents, SDC and MDC). Subsequent to 
compiling a list of plant species and environmental variables, we used various ordination methods 
to assess the relationships between vegetation coverage and environmental factors. Detrended 
correspondence analysis (DCA) was first utilized for ordination and the length of gradients was 
measured. Next, redundancy analysis (RDA) was selected due to the length of gradients (4.20.4; 
Table 1) (Ter Braak, 1985). The species responses against edaphic factors were assessed using 
plant species response curves derived from the GLM method (Salama et al, 2016). The 
relationships between vegetation coverage and edaphic factors were explored using these curves. 
The CANOCO v4.5 software was employed to perform ordination and draw the plant species 
response curves in relation to edaphic factors (Magurran, 2003). Moreover, cluster analysis of 
vegetation types in response to edaphic factors was performed using the Sorenson similarity 
indices and complete-linkage method embedded in the PAST v2.12 software. The stepwise 
multiple regression method was also used to explore the possible correspondence between SDC 
and MDC, as indices of soil aggregate instability and soil properties. 


Table1 Results of DCA (detrended correspondence analysis) of vegetation coverage and edaphic factors in the 
Great Salt Desert, central Iran 


Axis 1 2 3 4 
Eigenvalue 1.00 0.90 0.46 0.27 
Length of gradient 0.40 2.52 1.69 2.21 
Cumulative variance of plant species (%) 29.70 56.60 70.50 78.60 
Cumulative variance of soils (%) 33.70 54.80 0.00 0.00 


3 Results 


3.1 Density and coverage of vegetation communities 


The highest and lowest values of vegetation coverage were recorded for T. aphylla (68.196) and 
Tamarix aphylla-Halocnemum strobilaceum (2.6%) communities, respectively. Density of T. 
aphylla community ranged from high along t the streams to sparse with a low vegetation coverage. 
The highest plant density (33,000 stands/hm?) was observed in Phragmites australis-Halocnemum 
strobilaceum community and the least plant density (700 stands/hm?) belonged to Tamarix 
aphylla-Seidlitzia rosmarinus community (Table 2). 


3.2 Relationships between vegetation coverage and edaphic factors 


Soil physical and chemical properties measured in vegetation communities are presented in Table 
3. The highest EC value of 179.1 dS/m and pH value of 8.9 were observed in the S5 (T. aphylla) 


Table2 Elevation, vegetation coverage and plant density of different vegetation communities in the Great Salt 
Desert, central Iran 


Elevation Vegetation Plant density 


Number Vegetation community Site m) Oei (stands/hm?) 
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(%) 
: ; ; Hoz-e 
Si Phragmites australis-Halocnemum strobilaceum Soltan 818 6.5 33,000 
So Halocnemum strobilaceum-Tamarix aphylla Moreh 822 10.8 2300 
S3 Alhaji mannifera-Seidlitzia rosmarinus Haze 822 61.1 21,000 
Soltan 
S4 Tamarix aphylla-Halocnemum strobilaceum Moreh 824 2.6 3000 
Ss Tamarix aphylla Masileh 831 5.8 2000 
S6 Artemisia sieberi-Stipagrostis plumosa Moreh 850 2.8 10,000 
S; Artemisia sieberi Hees 855 3.4 6000 
Soltan 
; Hoz-e 
Ss Halostachys caspica-Halocnemum strobilaceum Soltan 817 11.5 4000 
So Halocnemum strobilaceum Hoze 819 30.7 10,000 
Soltan 
Sio Halocnemum strobilaceum-Seidlitzia Hoz-e 820 32.9 8600 
rosmarinus Soltan 
Su Halocnemum strobilaceum-Tamarix aphylla Moreh 800 3.4 2100 
; Hoz-e 
Si Tamarix aphylla Soltan 847 68.1 2900 
SB Tamarix aphylla-Seidlitzia rosmarinus Masileh 840 7.6 700 
Su Artemisia sieberi Masileh 855 4,7 9000 
Sis Halocnemum strobilaceum Masileh 819 6.9 3000 


and S,5 (H. strobilaceum) communities, respectively. The lowest SDC percentage was observed 
in Alhaji mannifera-Seidlitzia rosmarinus community (8.8%), highlighting the most stable soil 
structure with high sand but low clay contents in this region. MDC percentage of Artemisia 
sieberi-Stipagrostis plumosa community was as low as 11.1%. 

The ordination results obtained by RDA method showed that the study region encompassed 
three separate communities (Fig. 2). The first community included A. sieberi, S. plumosa and P. 
harmala, which were influenced by the soil sand, gypsum and carbonate contents. The second 
community consisted of H. strobilaceum, S. rosmarinus, H. caspica, P. australis and A. 
mannifera, whose presence was associated with high amounts of SDC and soil silt content. In the 
third community, the presence of T. aphylla was attributed to high levels of clay, soluble Na* 
concentration, SAR and EC. T. aphylla community was negatively correlated with soil sand 
content and the increment in soil soluble Na* concentration led to an increase in SAR. Moreover, 
EC, soluble Mg” concentration, SAR and clay content were relatively high in this vegetation type 
and showed positive correlations with the presence of T. aphylla. A. sieberi community was 
highly influenced by the percentages of sand and gypsum in the soils. H. caspica, H. strobilaceum 
and S. rosmarinus were positively linked to SDC and MDC and negatively to CCE content. P. 
australis and A. mannifera had positive correlations to soil silt content. Communities dominated 
by species such as S. plumosa and P. harmala showed positive correlations with CCE content 
(Fig. 2). Soluble Na* concentration, EC, CCE, SDC, clay, silt and sand contents were found to be 
the most important edaphic factors discriminating vegetation communities (Fig. 2). 

SDC and MDC were significantly and positively correlated with each other (r=0.61), while 
negatively with sand content with r values of —0.54 and —0.52, respectively (Table 4). The 
correlation analysis also revealed strong negative relationships between CCE content and 
dispersible clay contents (SDC and MDC) with r values of —0.70 and —0.84, respectively (Table 
4). SAR showed positive correlations with soluble Na* concentration (770.65) and EC (r=0.50). A 
negative relationship was found between EC and soil sand content (r= —0.42) (Table 4). 


MDC 
(%) 
23.2 
15.4 
26.5 
34.3 
26.0 
11.1 
42.2 
33.2 
50.8 
46.4 
67.1 
38.6 
25.0 
28.7 
22.6 


8.8 
34.3 
25.9 
11.1 
42.2 
33.2 
43.2 
46.4 
60.0 
30.8 
25.0 
13.3 
14.9 
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20.3 
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Fig. 2 Redundancy analysis (RDA) orientation diagram of different vegetation communities with environmental 
variables. EC, electrical conductivity; SAR, sodium adsorption ratio; OM, organic matter; SDC, spontaneously- 
dispersible clay; MDC, mechanically-dispersible clay; CCE, calcium carbonate equivalent. Mg”*, magmatism; 
Na*, sodium; Ca, calcium. S;—S15 represent the fifteen vegetation communities. 


Table 4 Correlation coefficients between edaphic factors through RDA 


Index Sand — Sit Clay OM CCE Gyp EC pH SAR Na Ca^ Mg* SDC MDC 
Sand — 1.00 


Sit — -0.88 . 1.00 

Clay  -030' 0.00 1.00 

OM 002 0.04 0.00 1.00 

CCE  -020 0.20 038  -036 1.00 

Gyp 028  -015  -030  -0.08 0.14 1.00 

EC  -042' 0.40 0.10 -0.30 -0.50° -0.30 1.00 

pH = -037  -0.54 J 030 0.10 025  -0.80 0.60° 1.00 

SAR  -023  -0.16 0.67  -023  -027  -030 0.50 022 1.00 
Nat  -040 0.29 0.26  -032  -040' -024 0.90 0.43"  0.65' 1.00 

Ca* -0.20 0.51 -0.60° -0.17  -038' -0.06 0.60 (0.90  -029 043° 1.00 

Mg* -041" O61 ' -030  -013 -0.50 -0.10 0.70 0.88  -022 0.62 078 ' 1.00 


0.22 0.40" 019 1.00 
0.33! 0.03 020 0.61” 1.00 


SDC -0.54" 001 -040° -0.07 -0.70° 0.57 019 0.54" 0.50 


* 


MDC -0.52" 0.11  -010 0.06  -0.84  -0.64 -0.25 031° 0.50 


Note: OM, organic matter; CCE, calcium carbonate equivalent; Gyp, gypsum; EC, electrical conductivity; SAR, sodium adsorption ratio; 
Na", sodium; Ca^, calcium; Mg”*, magmatism; SDC, spontaneously-dispersible clay; MDC, mechanically-dispersible clay. ', P<0.05 
level. 

The multiple regression model between dispersible clay content and soil properties revealed 
that SDC negatively linked with CCE (P«0.05, R°=0.477; Fig. 3a). MDC also had a negative 
relationship with CCE (P<0.05, R°=0.334; Fig. 3b). Results highlighted the importance of CCE as 
a predictor of SDC and MDC in the region. 

Cluster analysis based on vegetation coverage obtained three distinctive vegetation 
communities at the similarity level of 60% (Fig. 4). The first community belonged to carbonate- 
tolerant and coarse-textured soils-like species (psammophytes), which are found in desert areas 
with high soil sand contents. Species in the second community were found in areas with high soil 
silt contents. T. aphylla, as the most salinity- and sodicity-tolerant species, was categorized as the 
only member of the third community that grew in highly saline soils (Fig. 4). 


i; 70 ^ 
60- * SDC-68.19-1.935CCE (8) e MDC=64.75-1.546CCE (b) 
50}  R°=0.477, P<0.05 60- — q-0,334, P<0.05 
& 404 & ° 
Q Q 404 
A 304 Q 
n Z 30. 
204 3] 
e 
107 T T T T 104 ° 
I5 20 25 30 15 20 25 30 
CCE (%) CCE (%) 


Fig. 3 Relationships of spontaneously-dispersible clay (SDC, a) content and mechanically-dispersible clay 
(MDC, b) content with calcium carbonate equivalent (CCE) 
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Fig. 4 Cluster analysis based on vegetation coverage in the Great Salt Desert, central Iran 


3.3 Plant species response curves against edaphic factors 


Plant species differently respond to edaphic factors. The coverages of S. rosmarinus, T. aphylla 
and H. strobilaceum showed an upward trend with increasing gypsum content up to 1046 and then 
a downward one with gypsum content greater than 1546 (Fig. 5a). On the contrary, coverage of A. 
sieberi decreased with increasing gypsum content up to 10% and then started to increase with 
gypsum content greater than 15% (Fig. 5a). Coverage of T. aphylla increased while that of A. 
sieberi increased with increasing EC up to 150 dS/m. Coverage of H. strobilaceum and S. 
rosmarinus increased with increasing soil EC and started to decrease when EC values was greater 
than 100 dS/m (Fig. 5b). Coverages of S. rosmarinus, T. aphylla and H. strobilaceum increased 
with increasing soil sand content up to 50% and reduced when sand content was greater than 
60%. Coverage of A. sieberi showed a constant trend with increasing sand content and then had 
an upward trend when sand content was above 6096 (Fig. 5c). With increasing soil clay content, 
coverages of A. sieberi and H. strobilaceum decreased, whereas T. aphylla showed an upwardly 
trend. S. rosmarinus exhibited no discernible trend in relation to soil clay content (Fig. 5d). 


4 Discussion 


This study explored the relationships between vegetation coverage and soil characteristics in the 
Great Salt Desert, central Iran. The composition and distribution of vegetation communities in the 
study area were highly related to the ability of individual species to cope with flooding, salinity 
and sodicity, and nutrient limitation. For example, we found that the distribution of mono-specific 
communities (e.g., T. aphylla and H. strobilaceum) is related to specific soil physical and 
chemical gradients in the central Iranian deserts. Apaydin et al. (2009) also reported the 
occurrence of stable species around salt marsh of Black Sea. Significant zonation of saltmarsh 
plants appears in the Hoz-e Soltan Playa due to sharp physical gradients and species tolerance to 
harsh physical conditions (Sahragard and Chahouki, 2015). Edaphic factors and flooding are often 
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Fig. 5 Plant species response curves from generalized linear model (GLM) relating coverage of vegetation 
communities to edaphic factors. Ta.ap, Tamarix aphylla, Ha.st, Halocnemum strobilacum; Se.ro, Seidlitzia 
rosmarinus; Ar.si, Artemisia sieberi. 


reported to determine the establishment and zonal patterns of plant species around salt wetlands 
(Apaydin et al., 2009). This study found significant relationships between composition and 
distribution of vegetation communities and a series of edaphic factors including EC, soluble Na* 
concentration, SAR, CCE, sand and silt contents in the Great Salt Desert, central Iran , which is 
consistent with the result of Quevedo and Frances (2008), who assessed the relationship between 
soil and vegetation coverage in an arid and semi-arid region and found that soil texture was the 
most important contributing factor to establishment of plant species. The plant species response 
curves to environmental gradients in the study area provide a means for identifying the species 
fundamental and realized niches (Austin, 2005). It should be noted that the abundance of 
dominant plants in saline landscapes of central Iranian deserts is partly explained by the 
environmental gradients summarized in this study. Competition or facilitation between species 
and random recruitment should be explored to give a more accurate interpretation of vegetation 
changes as suggested by Zhang and Wang (2016). 

Hoz-e Soltan Playa, Moreh Plain and Namak Lake are among the lowest points in the center of 
the Iranian deserts, and facing to inlet and outlet discharged water in various seasons. The 
temperature fluctuations and flooding-drying cycles increase soil EC and soluble Ca”, Na* and 
Mg" concentrations in some vegetation communities in these areas. Soluble ions concentrations 
increase during summer as compared with spring and winter and then remarkably decrease during 
autumn season (Apaydin et al., 2009). Higher EC values in H. strobilaceum-T. aphylla, T. 
aphylla-H. strobilacum, T. aphylla and H. strobilaceum communities can also be related to the 
underground movement of saline waters from the Namak Lake towards the Masileh Plain and the 
subsequent upward transfer of salts to the soil surface due to evaporation and capillary rise. 
Furthermore, dam constructions in the upper zone of this area caused a sharp decrease in water 
input to this area compared with the past. In warm seasons, soil surface has undergone a 
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capillarity phenomenon due to the high temperature, causing increases in carbonate content, 
gypsum content, Na* and Mg” concentrations (Asrari et al., 2012). 

Relatively strong negative correlations were observed between dispersible clay content and 
CCE in the study area. These relations can be interpreted as follows. Carbonates are cementing 
agents for soil aggregation, and they release divalent cations (i.e., Ca^* and Mg”) in the soil 
solution that can negatively flocculate charged clay particles by cationic bridges. Therefore, 
increasing CCE leads to a decrease in aggregate instability as characterized by the dispersible clay 
content (Fig. 3). It should be noted that in general, soil salinity improves soil aggregation by 
contracting the diffuse double layer surrounding clay particles. In our study, however, the 
concentration of soluble Na* was higher than those of Mg” and Ca”, and therefore, the sodicity 
could minimize the effect of soil EC, leading to the disintegration of soil aggregates and 
increasing soil dispersible clay content, which is in agreement with the results of 
Gonzalez-Alcaraz et al. (2014) and Yang et al. (2016). 

Plant species with potential to tolerate the stress of drought environment can survive in these 
extreme environments (Kent, 2011). Hence, it is highly important to protect soil against erosion 
by identifying resistant species to such environment. According to the results obtained from the 
plant species response curves and RDA analysis, T. aphylla is the most resistant plant to 
unfavorable soil conditions in the study region. Moreover, A. sieberi can thrive in the region with 
high levels of carbonates. Plant species such as S. rosmarinusm and H. caspica can tolerate only a 
narrow range of environmental conditions and disappear beyond their edaphic thresholds 
(Hejcmanova-Nezerkova and Hejcman, 2006; Kent, 2011). The structural and functional 
responses of plant species to various environmental conditions can be well assessed using an 
appropriate multivariate analysis approach (Kleyer et al., 2012). The results of this study may 
lead to new approaches for more effective solutions to the saline environment of central Iranian 
deserts. 


5 Conclusions 


Our results found that vegetation communities in the study area were classified into three 
distinctive vegetation communities, which included carbonate-tolerant and coarse-textured- 
dependent species, high soil silt and organic matter content-dependent species, and salinity- and 
sodicity-tolerant species. We also highlighted the importance of soil carbonate content as a 
predictor of dispersible clay content in these arid and saline areas. 

Due to the anthropogenic effects and frequent drought events, central Iranian deserts are often 
experiencing severe degradation in soil and vegetation. In order to avoid permanent loss of the 
plant diversity, it is necessary to implement conservation strategies in these fragile ecosystems. 
Thus, our data provide a valuable baseline for future monitoring studies regarding the possible 
changes in plant species distributions, and can be used to facilitate reclamation and restoration 
projects. 
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